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摘 要 : 地 球 辐射 带 中 的 高 能 带电 粒子 是 引起 航天 器 材料 和 器 件 性 能 退化 甚至 失效 
的 主要 空间 环境 因素 . 因此 ,航天 器 设计 中 所 采用 的 辐射 带 模型 的 准确 程度 对 于 航天 器 的 生存 
能 力 和 航天 任务 的 完成 质量 至 关 重 要 . 在 利用 我 国 自 主 辐射 带 高 能 粒子 探测 数据 进行 的 辐射 
带 建 模 中 , 探测 数据 的 空间 网 格 化 是 一 项 非常 重要 的 工作 . 介绍 了 我 国 辐射 带 探测 数据 的 情 
况 ,以 及 辐射 带 建 模 的 方法 和 步骤 ; 重点 研究 了 不 同 插值 方法 在 低地 球 轨道 ( LEO) 空间 辐射 带 
建 模 数据 网 格 化 中 的 应 用 ,并 开展 了 误差 分 析 . 研究 结果 表明 : 在 各 种 常用 的 插值 方法 中 , 反 距 
离 加 权 法 、 自 然 邻 点 法 和 最 近邻 点 法 适合 工程 化 应 用 . 其 中 , 反 距 离 加 权 法 生成 的 数据 网 格 对 
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粒子 通 量 的 反 演 结果 精度 最 高 , 该 方法 采用 低 阶 距离 时 得 到 的 反 演 结果 更 为 合理 . 
X WE 词 : 空间 环境 ; 辐射 带 ， 数据 网 格 ; 插值 ， 反 距离 加 权 法 ; 自然 邻 点 法 ; 最 近邻 
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地 球 外 层 空间 存在 着 一 个 区 域 ,其 中 充满 地 
磁场 捕获 的 高 能 质子 和 电子 , 这 个 区 域 被 称 为 地 
球 辐射 带 ( 以 下 简称 为 辐射 带 ) " . 辐射 带 中 的 质 
子 和 电子 能 量 较 高 ” ,能 够 引起 航天 器 材料 和 器 
件 性 能 退化 甚至 失效 . 在 地 球 空间 运行 的 绝 
大 多 数 航天 器 都 要 或 多 或 少 地 穿越 辐射 带 , 遭遇 
高 能 粒子 辐射 . 因此 ,在 航天 任务 的 设计 过 程 中 ， 
准确 地 定量 估算 航天 器 所 处 的 空间 辐射 环境 对 于 
航天 器 设计 的 优化 、 运 行 安 全 的 保障 和 任务 的 完 
成 都 至 关 重 要 . 

目前 国际 和 国内 航天 界 普 吾 采 用 NASA 编制 
的 质子 模型 AP8 和 电子 模型 AE8 来 计算 辐射 带 
质子 和 电子 的 通 量 , 并 以 此 为 基础 评估 辐射 带 对 
航天 器 的 影响 ”< . 目前 我 国 可 获得 的 AP8/AE8 
模型 的 最 新 版 本 分 别 完成 于 1976 年 和 1983 ££, 
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这 两 个 版 本 使 用 的 都 是 20 世纪 70 年 代 以 前 的 数 
2" ,观测 资料 比较 陈旧 ,模型 计算 值 与 实际 观 
测 值 经 常 存在 较 大 偏差 ”3 ,使 得 这 两 个 版 本 
的 模型 应 用 到 当前 航天 任务 中 存在 时 效 局 限 性 . 
NASA 在 后 期 对 APS/AES 模型 进行 了 多 次 修订 
和 补充 ,形成 了 多 种 修正 版 本 ,并 且 于 2006 年 联 
合 美国 多 家 单位 开始 开发 用 于 替代 APS/AES 的 
辐射 带 质 子 和 电子 的 新 模型 : APO /AEO " ,但 由 
于 技术 封锁 ,我 国 目 前 还 无 法 获得 APS / AES. WE 
正版 本 和 AP9/AF9. 因此 ,我 国 自主 研发 辐射 带 
模型 对 于 我 国航 天 工程 的 可 持续 发 展 具 有 现实 和 
战略 的 意义 . 地 面 高 度 为 300 ~ 1000 km 的 轨道 空 
间 即 低地 球 轨道 ( Low-Earth Orbit, LEO) 空间 是 大 
量 应 用 卫星 运行 的 区 域 .因此 ,这 一 区 域 是 辐 
射 带 建 模 中 最 重要 的 区 域 之 一 . 
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1 我 国 辐射 带 监测 数据 情 ; 

我 国 从 20 世纪 80 年 代 开始 对 空间 环境 展开 
了 大 规模 探测 ,其 中 LEO 空间 辐射 带 高 能 粒子 分 
布 是 探测 的 重点 ”. 截至 到 目前 为 止 ,已 发 射 卫 
星 中 搭载 有 高 能 粒子 探测 器 的 卫星 共有 50 余 颗 ， 
其 中 LEO 卫星 超过 30 颗 . 部 分 高 能 粒子 探测 仪 
器 的 卫星 轨道 .观测 对 象 . 能 量 范围 .通道 数 和 观 
测 时 段 如 表 1 所 示 . 这 些 仪器 的 设计 原理 相同 ,者 


常 峥 ,等 : 基于 数据 网 格 化 方法 的 低 轨 辐射 带 建 模 技术 2289 


为 半导体 型 探测 器 . 

不 同 仪器 获得 的 探测 数据 在 联合 使 用 前 , 必 
须 经 过 交叉 定 标 “ . 因此 在 建 模 过 程 中 , 研究 人 
员 对 表 1 所 列 数据 进行 了 同化 处 理 , 并 以 此 为 基 
人 础 建立 了 LEO, 中 地 球 轨道 ( Medium 下 arth Orbit, 
MEO) 和 同步 轨道 ( Geosynchronous Orbit, GEO) pr 
射 环境 数据 库 , 数据 库 包含 了 超过 一 个 完整 太阳 
活动 周 ( 11 a) 的 LEO 空间 辐射 带 数 据 , 可 以 用 于 
LEO 空间 辐射 带 建 模 . 


表 1 高 能 粒子 探测 仪器 参数 
Table 1 High energy charged particle detector parameters 


卫星 名 称 轨道 观测 对 象 ” 能 量 范围 /MeV 通道 数 观测 时 段 
zt 质子 能 谱 2.9 ~300 10 
和 四 号 、 号 近 圆 轨道 ,高 度 340 km, 倾角 42° A à 200242—2004-05 
神舟 四 号 HAAS 近 圆 轨道 ,高 度 m, 倾角 电子 能 谱 ar 9 
a "— E 质子 能 谱 >1.5 8 
天 宫 一 号 近 圆 轨道 ,高 度 380 km, 倾角 43 ms P $ 201140 至 今 
= 质子 能 谱 3 ~300 6 
风云 三 号 A 星 、.B 星 、C 近 圆 轨道 ,高 度 830 km, tif 98° 2008-05 ES 
风云 三 号 A 星 .B 星 、C 星 近 圆 轨道 , 高度 my 倾角 9 电子 能 谱 015-57 $ 至 今 
= 、 ARE — 2.9 ~300 5 
风云 一 号 C 星 、.D 星 近 圆 轨道 , SE 860 km, 倾角 98° 质子 能 谱 1995-05 一 2010-06 
电子 通 量 >1.6 1 
ae 质子 能 谱 2.9 ~300 5 
x a 近 圆 轨道 ,高 度 860 km , 倾 鱼 98。 1999-95—2000-02 
实践 五 号 近 圆 轨道 ,高 度 my 倾角 9 电子 能 谱 015-57 999 
E 大 椭圆 轨道 ,近地点 200 km, 质子 能 谱 4 ~300 
R 号 一 1994-02—1994-08 
实践 四 号 远地点 36000 km, 倾角 28. 5? 电子 能 谱 0.5 ~4 5 a 2 
质子 能 谱 ”2.9 ~300 1997-06 一 1998-04 
风云 二 号 A 星 、B 星 静止 轨道 质子 通 量 >1.1 1 se de DOT 
电子 通 量 >1.4 1 
质子 能 谱 10 ~300 3 
风云 二 号 C 星 .D 星 工 星 静止 轨道 >0.35 200440 至 今 
= 4 
电子 通 量 d 
质子 能 谱 >4 6 
风云 二 号 下 道 2012-01 至 今 
风云 二 号 下 星 静止 轨道 mun oe i BS 


2 辐射 带 建 模 的 方法 

本 文 利 用 我 国 辐射 带 高 能 粒子 的 探测 数据 建 
立 的 自主 辐射 带 模型 和 NASA 的 AP8/AF8 均 为 
经 验 模 型 , 即 模型 的 基础 是 具有 特定 组 织 结 构 的 
数据 ,文中 将 这 种 具有 特定 组 织 结构 的 数据 定义 
为 表单 . 在 辐射 带 模 型 中 , 表单 指 的 是 能 描述 辐射 
环境 某 一 给 定时 刻 ` 给 定 空 间 环境 状态 下 在 所 有 
空间 维度 上 分 布 特征 的 粒子 通 量 数据 列表 . 
AP8/AE8 仅 含有 一 组 代表 太阳 活动 高 年 和 低 年 时 
对 应 高 能 粒子 通 量 平均 分 布 的 表单 ,不 能 给 出 辐 
射 带 随地 磁场 长 期 变化 而 引起 的 缓 变 , 也 不 能 
映 太 阳 活 动 周 内 的 各 种 周期 性 变化 和 各 种 突 发 性 
扰动 变化 吗 ,因此 称 其 为 静态 模型 ; 代替 它 的 
AP9 /AE9 只 能 通过 输出 结果 中 的 概率 间接 地 反 
映 各 种 变化 对 辐射 带 的 影响 呈 ,说 明 APO/AEO 
也 不 是 真正 意义 上 的 动态 模型 . 


本 文 试 图 建立 的 自主 辐射 带 模型 的 建 模 主要 
步骤 如 图 1 所 示 . 从 图 1 中 可 以 看 出 ,模型 输出 结 
果 是 通过 对 模型 数据 中 的 4 组 表单 进行 查询 、 计 
算 并 将 结果 进行 于 加 后 得 到 的 ,因此 这 4 组 表单 
是 模型 的 核心 内 容 . 这 4 组 表单 分 别 给 出 辐射 带 
的 宁静 状态 .长 期 缓 变 状 态 ` 周 期 变化 状态 和 扰动 
变化 状态 ,因此 自主 辐射 带 模型 属于 动态 模型 ,这 
也 是 该 模型 相 比 于 AP8/AES8 的 最 大 改进 之 处 . K 
单 的 建立 流程 如 图 2 所 示 . 同化 后 的 数据 在 空间 


i tn, f= SF ky dpe BS 
H 静态 平均 结果 | 


长 期 组 变数 据 表单 ES UIT 
| [ 司 二 性 变化 数据 表单 上 | 同期 性 变化 结果 


| [ 突 发 性 变化 数据 表单 上 一 [ 突 发 性 变化 结果 | 一 


J 图 1 自主 辐射 带 模 型 框架 


Fig.1 Framework of own radiation belts model 
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企 空 间 中 均匀 分 布 后 ,再 进行 分 类 处 理 , 生 成 各 类 
数据 库 . Alt, 空间 网 格 化 是 建 模 中 的 重要 步骤 . 


建立 表单 
空间 、 粒 子 人 射 方向 、 静态 数据 
z 能 谱 、 粒 子 种 类 N 表单 
i D i 1 
as [n] i] i 地 磁场 长 期 变化 i KHZ uy 
数 门 | 化 ^ = a M 数据 表单 
据 m 2 TW 
Gi ! 太阳 活动 周 变化 、 MITES 
| 季节 变化 、27 天 周期 变化 上 -| e e 
地 方 时 变化 sa 
WASP. 分 级 p. 突 发 性 变化 
UT Lie ae 


2 自主 辐射 带 模型 表单 建立 流程 
Fig.2 Procedure of creating data table in 


own radiation belts model 


3 种 轨道 空间 中 ,LEO 空间 卫星 最 多 , H LEO 2E 
间 包 含 了 多 个 不 同 高 度 的 轨道 ,而 MEO( 特 指 高 
度 为 22 000 km、 倾 角 为 55° 的 卫星 轨道 ) 空间 和 
GEO 空间 均 仅 包含 单一 轨道 的 若干 颗 卫 星 . 
此 ,LEO 空间 的 数据 网 格 化 成 为 建 模 中 的 重点 . 


3 LEO 空间 数据 网 格 化 


3.1 空间 数据 网 格 

鉴于 LEO 空间 观测 数据 的 分 布 特点 ,为 了 得 
到 如 图 3 所 示 的 空间 数据 网 格 , 本 文采 取 了 以 下 
处 理 方法 进行 数据 网 格 化 : 

1) 利用 11 个 轨道 高 度 上 的 卫星 观测 资料 ， 
完成 对 应 11 个 空间 球面 上 的 数据 网 格 化 . 对 于 同 
系列 卫星 观测 数据 的 使 用 原则 是 ,利用 其 中 部 分 
卫星 的 观测 资料 完成 插值 ,剩余 卫星 的 观测 资料 
用 于 对 插值 结果 进行 验证 . 

2) 11 个 空间 球面 以 外 的 空间 ,利用 11 个 空 
间 球 面 的 数据 网 格 ,结合 大 椭圆 轨道 卫星 的 观测 
资料 通过 插值 获得 ; 11 个 空间 球面 之 间 的 区 域 采 
用 内 插 , 高 度 范围 为 300 ~ 340 km、860 ~ 1000 km 
的 区 域 采用 外 插 . 


3 LEO 空间 数据 网 格 
Fig.3 Data grid of LEO space 
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3.2 可 选用 的 插值 方法 

空间 网 格 点 上 的 粒子 通 量 大 小 是 通过 对 观测 
数据 插值 获得 的 , 而 插值 中 最 关键 的 步骤 是 利用 
同一 空间 球面 上 观测 数据 进行 插值 得 到 该 球面 上 
的 网 格 点 通 量 . 常用 的 插值 方法 包括 反 距 离 加 权 
法 ( Inverse Distance Weighting, IDW) i BK 
法 ( Kriging ) 09. 自然 邻 点 法 ( natural neigh— 
bor) pa \、 最 近邻 点 法 ( nearest neighbor) ol 多 元 
回归 法 ( polynomial regression) 中、 最 小 曲率 法 
( minimum curvature) 5324 C] 连续 5 次 双 变 量 插 
值 法 ( quintic) “”、 改 进 谢 别 德 法 ( modified 
Shepard’ s) 0723 e S EEG ( Radial Basis Func- 
tion, RBF) pal 、 三 角形 剖 分 法 ( triangulation) ds 
等 . 这 些 方法 在 基本 原理 、 数 学 模型 .计算 效率 等 方 
面 各 不 相同 , 但 在 使 用 中 一 般 都 需要 输入 辅助 参 
数 , 例 如 : 克 里 格 法 需要 输入 所 采用 的 变异 函数 的 
XA DE REE REA” 改进 谢 别 德 法 需 
输入 节 画 数 的 两 个 影响 半径 ”; OR E] EE ERG SS 
要 确定 插值 使 用 的 径 向 基 瑟 数 ”. 

在 应 用 插值 方法 时 , 由 于 网 格 点 和 观测 数据 
分 布 在 球面 上 ,因此 插值 形式 应 为 球面 插值 ; 同时 
由 于 辐射 带 内 的 粒子 通 量 分 布 服从 宕 律 谱 只 ， 
此 计算 中 应 采用 对 数 插值 . 

使 用 空间 球面 上 的 网 格 点 通 量 可 以 反 演 出 该 
球面 上 任意 位 置 的 通 量 . 通过 对 反 演 结果 的 误差 
分 析 可 以 确定 网 格 点 通 量 的 精度 , 进而 确定 选取 
的 插值 方法 是 否 合 适 . 对 反 演 结果 的 误差 分 析 主 
要 是 对 反 演 值 ( P) 和 观测 值 ( 0) Z#(D =P - Q) 
进行 分 析 . 误差 分 析 包 括 两 部 分 : 计算 各 类 误差 成 
分 在 总 误差 中 的 比例 和 计算 误差 大 小 外 . 本 文采 
用 均 方 误差 ( Mean Square Error, MSE) 来 描述 误差 
比例 . MSE 包括 系统 均 方 误差 ( systematic MSE) 
MSE, F0 3E 3& £t 15 75 1x Æ ( unsystematic MSE) 
MSE, ,其 中 MSE, 由 3 部 分 组 成 ,包括 MSE 中 的 
偏 移 成 分 ( additive MSE) MSE, 、 缩 放 成 分 ( propor- 
tional MSE) MSE, 和 MSE,. MSE, 互相 依赖 成 分 
( interdependence MSE) MSE,. 以 上 误差 的 计算 方 
法 由 式 (1) ~ 式 (6) 给 出 : 


3 UR -0)' 
_ i=l 
MSE = N (1) 
MSE, =a’ (2) 
N 
0; 
MSE, 2 (b - 1) = 
SE, - ( p (3) 
MSE, =2a(b - 1) O (4) 
MSE, = MSE, + MSE, + MSE, (5) 


第 12 其 

MSE, = MSE - MSE, (6) 

TAN 为 观测 点 数 ; P, 和 0; 分 别 为 第 i 点 的 计 

算 值 和 观测 值 ; 0 为 观测 值 的 平均 值 ; Ab 分别 

为 以 0 为 自 变 量 \ 以 P 为 因 变 量 采用 最 小 二 乘法 
拟 合 出 的 直线 的 截 距 和 和 斜率 . 

本 文采 用 RMSE、RMSE.、RMSE,、RMSE,、 
RMSE，、RMSE; 来 描述 误差 大 小 ,这 6 个 参数 值 的 
量 纲 与 粒子 通 量 相同 ,大 小 分 别 等 于 MSE.MSE, . 
MSE, . MSE, . MSE, . MSE, 的 平方 根 ,符号 与 MSE、 
MSE.、MSE,、MSE,、MSE,、MSE, 相同 . 同时 , 还 采 
用 了 一 致 性 系数 d 来 描述 P 和 0 的 一 致 程度 ,d 
的 计算 方法 为 


Y (P, a 0) 
d=1-—> (7) 
2, (IP: l+ lo; 
式 中 : Ps A 071 分别 由 式 (8) 和 式 (9) 给 
P;=P,-0 (8) 
0;20,-0 (9 
d 的 取 值 范围 在 0.0 ~ 1.0 之 间 ,d 越 接近 1.0 


表示 计算 值 和 观测 值 的 一 致 性 越 高 . 
3.3 插值 方法 的 应 用 

本 文通 过 对 轨道 高 度 为 630 km 的 两 颗 卫 星 
的 观测 数据 的 处 理 来 说 明 空 间 球面 上 的 数据 网 格 
化 方法 及 精度 评估 标准 . 这 两 颗 卫星 为 同系 列 卫 
星 , 各 自 搭 载 了 一 台 由 中 国 科 学 院 空 间 中 心 空 间 
环境 探测 研究 室 研 制 的 高 能 粒子 探测 器 ,两 台 探 
测 器 的 工作 原理 、 制 造 工艺 和 技术 指标 完全 相同 . 
该 探测 器 可 以 对 6 个 能 道 的 高 能 质子 通 量 和 1 个 
能 道 的 高 能 电子 通 量 同 时 进行 观测 . 选 定 的 观测 
数据 为 2010 年 8 月 29 日 至 10 月 10 日 期 间 能 量 
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范围 为 3 ~5 MeV 的 质子 通 量 , 该 时 间 段 内 空间 环 
境 始终 保持 宁静 , 得 到 的 观测 数据 中 不 包含 突 发 
性 空间 环境 扰动 引起 的 质子 通 量变 化 . 

数据 处 理 的 具体 步骤 如 下 : 

1) 利用 其 中 一 颗 卫星 的 数据 ,选取 一 种 插值 
方法 ,给 定 该 方法 所 需 的 各 项 参数 值 ,得 到 卫星 所 
处 空间 球面 上 的 网 格 点 通 量 , 网 格 划 分 以 经 度 、 纬 
度 为 单位 ,每 1° 作 为 一 个 网 格 . 

2) 使 用 网 格 点 通 量 反 演 另 一 颗 卫 星 各 观测 
位 置 的 通 量 , 反 演 采用 趋势 面 插值 法 , 即 在 网 格 点 
中 选取 3 个 点 ,此 3 个 点 形成 的 三 角形 是 可 以 包 
围 采样 位 置 的 最 小 三 角形 ,将 网 格 点 的 经 度 、 纬 度 
作为 x、y 坐标 , 通 量 作为 z 坐标 , 拟 合 出 一 个 趋势 
面 方程 ,利用 该 方程 计算 出 观测 位 置 的 通 量 . 

3) 将 反 演 值 与 观测 值 进行 比 对 ,得 到 所 选 定 
插值 方法 在 采用 选 定 参 数 时 的 误差 分 析 结 果 . 

在 第 3.2 节 给 出 的 各 种 插值 方法 中 , 反 距 离 
加 权 法 、 自 然 邻 点 法 、 最 近邻 点 法 、 克 里 格 法 和 径 
向 基 事 数 法 支持 球面 插值 , 克 里 格 法 和 笃 向 基 事 
数 法 在 使 用 中 需要 较 多 的 人 工 干预 才能 得 到 精度 
较 高 的 结果 ,不 适 于 工程 应 用 . 因此 ,在 实际 建 模 
中 本 文选 择 了 反 距离 加 权 法 、 自 然 邻 点 法 和 最 近 
邻 点 法 进行 对 比 . 表 2 给 出 了 这 3 种 插值 方法 处 
理 结果 的 误差 分 析 对 上 比 ,改变 反 距离 加 权 法 的 输 
入 参数 会 得 到 不 同 的 处 理 结果 , 表 2 给 出 的 是 其 
误差 分 析 最 优 的 结果 . 从 表 2 给 出 的 结果 可 以 看 
出 , 反 距 离 加 权 法 在 最 优 状 态 下 得 到 的 处 理 结 
的 精度 高 于 自然 邻 点 法 和 最 近邻 点 法 ,原因 在 于 
其 输入 参数 优化 处 理 能 够 大 幅 减 小 系统 偏差 
( RMSE,) 中 的 缩放 成 分 ( RMSE.) , 自然 邻 点 法 和 
最 近邻 点 法 没有 可 调整 的 输入 参数 , 反 距 离 加 权 
法 有 两 个 输入 参数 可 以 调整 ,其 计算 原理 为 


表 2 采用 反 距离 加 权 法 、 自 然 邻 点 法 和 最 近邻 点 法 处 理 结果 的 误差 分 析 对 上 比 
Table 2 Error analysis comparison of processing results of IDW, natural neighbor and nearest neighbor 
方法 RMSE/ RMSE, / RMSE, / 
( cm 2 £) (cm~? « sr“! * s7!) (cm~? * sr7! * s7!) 
自然 邻 点 法 -0.69511 1.25167 71.68575 46. 598 82 54.473 82 
最 近邻 点 法 -0.48722 1.25192 73. 888 11 49. 802 22 54. 58197 
反 距 离 加 权 法 -2.19672 1.01709 48. 647 80 48. 483 42 3.995 68 
方法 RMSE, / RMSE, / RMSE, / MSE, /MSE MSE MS 
(em? * s^! * s!) (em? * s^! + s7!) (em? *sr^! * s!) 
自然 邻 点 法 0.69511 54.58076 -3.48511 0.42256 0.57744 
最 近邻 点 法 0.48722 54. 657 84 -2.91982 0.45431 0. 545 69 
反 距 离 加 权 法 2. 19672 3.70774 -1.61477 0.99325 0.00675 
方法 MSE, /MSE MSE, /MSE MSE, /MSE d 
自然 邻 点 法 0.00009 0.57971 -0.00236 0.97811 
最 近邻 点 法 0. 00004 0.54721 -0.00156 0.97678 
反 距离 加 权 法 0.00204 0. 005 81 -0.001 10 0. 987 48 
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j= + 一 (10) 
Xv 

式 中 :$ 为 待 插 值 点 的 属性 值 , 即 网 格 点 通 量 ; 
v; 为 参与 插值 的 离散 点 属性 , 即 观测 位 置 的 通 量 ; 
d; 为 观测 位 置 和 网 格 点 的 距离 ,因此 , 可 以 通过 
调整 距离 的 阶 数 p 和 观测 数据 的 范围 ”获得 不 同 
的 插值 结果 ,而 在 实际 应 用 中 一 般 是 通过 改变 搜 
索 半径 7 来 控制 参与 插值 的 观测 数据 的 范围 . 

在 计算 中 ,将 > 的 范围 确定 为 1 ~22013. 14 km, 
其 中 1 km 为 观测 点 至 网 格 点 距离 的 最 小 值 ， 
22013. 14 km 为 卫星 轨道 球面 上 大 圆 的 半 周 长 . 在 
此 范围 下 , 本文 获 得 了 p 从 1 ~10 的 处 理 结 
图 4 和 图 5 所 示 分 别 为 不 同 p EX (& P BS d 和 
RMSE 的 对 比 . 从 图 4 可 以 看 出 , 当 r < 1000 km 
时 ,不 同 p 取 值 下 的 结果 基本 相同 ; 4 >1000km 
时 ,在 p =1 和 p=2 的 情况 下 d 会 随 着 7 的 增加 迅 
速 减 小 ,在 p 宇 3 的 情况 下 d 基本 保持 不 变 ,仍然 
维持 较 高 水 平 . 图 5 所 示 的 RMSE 的 变化 规律 与 
d 类 似 , 在 p=1、2 和 3 的 情况 下 最 小 RMSE 值 均 
小 于 p 宇 4, 但 获得 最 小 RMSE 的 > 各 不 相同 .图 4 
和 图 5 的 结果 说 明 , TE p 2 4 的 情况 下 , 当 7 > 
65 km 时 ,不 同 p f& F d 和 RMSE 的 变化 规律 几乎 
完全 相同 ,说 明 此 种 情况 下 单纯 增 大 7 不 能 影响 
插值 结果 的 精度 , 这 是 由 反 距 离 加 权 法 的 原理 决 
定 的 . 从 式 ( 10) 中 可 知 ,p 增 大 会 增加 距离 网 格 点 
较 近 的 实测 点 在 插值 结果 中 的 权重 . 因此 , 4 r BX 


10 000 


2015 年 
值 使 得 搜索 范围 包含 全 部 对 插值 结果 有 主要 作用 
的 实测 点 后 , 增 大 > 基本 不 会 改变 插值 结果 的 
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Fig.4 Comparison of processing results of d with 
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— p=l 

250 — p=2 
p=3 

200 — p^ 
— P= 


RMSE/(cm^'sr *s^') 


0 eee ae MERI. 
10 100 ] 000 
r/km 


5 反 距 离 加 权 法 取 不 同 p 时 处 理 结果 的 RMSE 对 比 
Fig.5 Comparison of processing results of RMSE with 
different p in IDW 


图 6 给 出 了 反 演 结果 和 观测 结果 的 对 比 , 反 
演 结果 采用 反 距 离 加 权 法 处 理 获 得 ,P 从 1 取 到 
10, 选 定 的 7 使 当前 p 下 反 演 结果 的 RMSE 最 小 . 
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6 ”质子 通 量 观测 结果 和 网 格 数据 反 演 结果 对 上 比 


Fig.6 Proton flux comparison of observation and inversion result by data grid 
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从 图 6 可 以 发 现 观测 结果 具有 以 下 缺陷 : 

1) 本 底 较 高 ,基本 维持 在 大 于 10/( cm * sr * s) 
的 水 平 ,这 是 由 探测 器 几何 因子 较 小 决定 的 . 

2) 通 量 范围 在 10 ~ 100/( em? * sr * s) 的 区 
间 内 ,观测 结果 抖动 比较 严重 ,这 是 受 仪器 本 底 噪 
声 和 观测 数据 采用 遥测 分 层 值 方式 采集 的 共同 影 
响 造 成 的 . 

理想 的 网 格 数据 给 出 的 反 演 结果 ,应 该 可 以 
对 观测 数据 中 出 现 的 以 上 缺陷 予以 修正 . 对 比 
6 给 出 的 各 阶 反 演 结果 ,可 以 发 现 它 们 有 具有 以 
下 特点 : 

1) 通 量 大 于 100/( em? * sr * s) 时 ,各 阶 反 演 
结果 与 观测 结果 一 致 性 均 较 高 . 

2) 通 量 小 于 100/( em? * sr * s) 时 , 低 阶 (p = 
l,p =2) 反 演 结果 具有 较 好 的 消除 抖动 的 效果 ， 
对 于 本 底 以 下 的 通 量具 有 较 好 的 外 推 效果 . 

因此 ,使 用 反 距 离 加 权 法 进行 处 理 , 当 距 离 阶 
EX p 选 定 为 1 或 2 时 ,得 到 的 反 演 结果 合理 程度 
更 高 ; 同时 , 根据 图 4 和 图 5 给 出 的 误差 分 析 结 
果 , 选 择 合适 的 搜索 半径 7, 可 以 使 反 演 结果 精度 
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本 文 说 明了 自主 辐射 带 模型 属于 动态 模型 的 
原因 ,指出 了 自主 模型 相 比 于 APS/AES 的 改进 之 
处 , 即 可 以 给 出 各 种 辐射 带 周 期 性 变化 和 扰动 变 
化 ; 讨论 了 辐射 带 自主 建 模 的 目标 和 思路 ,针对 
LEO 空间 数据 网 格 化 问题 ,根据 目前 我 国 自主 探 
测 数 据 的 实际 情况 ,使 用 了 空间 分 层 的 方式 对 每 
一 个 卫星 所 处 空间 球面 上 的 数据 分 别 进行 处 理 ; 
在 生成 数据 网 格 的 插值 计算 中 , 根据 各 常用 插值 
方法 的 特点 ,选择 反 距 离 加 权 法 、 自 然 邻 点 法 和 最 
近邻 点 法 作为 考察 的 插值 方法 ,通过 对 网 格 数据 
有 反 演 结果 的 误差 分 析 和 和 与 观测 结果 的 比 对 ,确定 
有 反 距 离 加 权 法 具有 更 高 的 精度 , 且 反 距离 加 权 法 
在 采用 低 阶 距离 的 情况 下 生成 的 数据 网 格 得 到 的 
反 演 结果 具有 更 高 的 合理 性 ,并 可 以 对 实测 数据 
中 的 一 些 缺 陷 进行 弥补 . 

本 文 讨 论 的 LEO 空间 数据 网 格 化 方法 的 评 
价 标准 建立 在 误差 分 析 的 基础 上 , 是 一 种 数值 分 
析 的 方法 . 实际 上 , 辐射 带 的 形成 有 着 复杂 的 物理 
机 制 , 其 内 部 高 能 带电 粒子 的 分 布 也 服从 相应 的 
物理 规律 ,因此 辐射 带 建 模 不 能 脱离 物理 理论 的 
指导 . 在 后 续 工作 中 需要 重点 改进 的 地 方 是 将 物 
理 理论 与 数值 分 析 方 法 相 结 合 ,利用 物理 理论 对 
数值 分 析 方 法 进行 完善 . 
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Modeling technology of radiation belt in LEO based on 
data gridding methods 
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Abstract: The energetic charged particle in the radiation belts of earth is the main environmental factor 
which can result in performance degradation and even failure of the material and device on spacecraft, so the 
accuracy of radiation belts model which is being used during design phase of spacecraft is very important to 
spacecraft performance of space mission and survivability. In the development of radiation belts model based 
on our own country’s data from exploration of radiation belts, the gridding of discrete data is a fundamental 
task. The key facts of data from Chinese exploration of radiation belts were introduced, the goal and road-map 
of modeling were described, the application of interpolation in gridding of modeling of radiation belts in low- 
earth orbit( LEO) was discussed, and the significant errors were analysed. The result of interpolation indicates 
that in mainly interpolation methods, inverse distance weighting ( IDW) , natural neighbor and nearest neigh- 
bor are most suitable for engineering calculation; the particle flux which is inverted from data grid produced by 
IDW has the highest precision, and the result from IDW will be more accurate when IDW uses small exponent 
of distance. 
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